Local-fields and disorder effects in free-standing and embedded Si nanocrystallites 
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The case study of a 32-atoms Si nanocrystallite (NC) embedded in a Si02 matrix, both crystalline 
and amorphous, or free-standing with different conditions of passivation and strain is analyzed 
through ab-initio approaches. The Si32/Si02 heterojunction shows a type I band offset highlighting 
a separation between the NC plus the interface and the matrix around. The consequence of this 
separation is the possibility to correctly reproduce the low energy electronic and optical properties of 
the composed system simply studying the suspended NC plus interface oxygens with the appropriate 
strain. Moreover, through the definition of an optical absorption threshold we found that, beside 
the quantum confinement trend, the amorphization introduces an additional redshift that increases 
with increasing NC size: i.e. the gap tends faster to the bulk limit. Finally, the important changes 
in the calculated DFT-RPA optical spectra upon inclusion of local fields point towards the need of 
a proper treatment of the optical response of the interface region. 



INTRODUCTION 



II. THEORETICAL METHOD 



C 

O 

o 



> 
-si- 
r- 
en 
p 

ON 

o 
o 



X 



Light emission in the visible spectral range and at room 
temperature has been observed in the last two decades 
in porous Si before and in nanostructured Si after. Si 
nanocrystallites embedded in a silicon dioxide matrix, 
in particular, have attracted much attention in recent 
years 1-3 . In these structures, Si acts as the low-energy 
gap material. The strong quantum confinement induced 
by the dielectric matrix enlarges the band-gap and in- 
creases the spatial overlap of the electron and hole wave 
functions, thus increasing the spontaneous emission rate 
and shifting the emission peak to energies higher than 
bulk Si band-gap. Despite the large number of exper- 
iments and theoretical studies, many question are still 
open concerning the nature of the optical properties of 
Si-NCs embedded in a silica matrix. Actually, too many 
parameters, like size, shape, strain and interface struc- 
ture, influence the spectral response of the system. The 
characterization of these quantities is thus of fundamen- 
tal importance in order to determine the best condition 
for light emission and optical gain. 

In this work, we present DFT calculations on a crys- 
talline and amorphous 32-atoms Si-NC, passivatcd with 
hydrogens, OH groups, and embedded in a Si02 matrix. 
The comparison between the embedded and suspended 
NCs can give insights on the role played by strain in 
the determination of the electronic and optical properties 
of the composed system, while the different passivation 
regimes can reveals the contribution given separately by 
the NC, the interface and the matrix to the properties of 
the whole system. 

The paper is organized as follows. A description of the 
theoretical methods and of the systems is given in section 
II. The electronic and optical properties of the considered 
systems are analyzed in section III while conclusions are 
presented in section IV. 



Ab-initio calculations of structural, electronic and opti- 
cal properties for a 32 atoms Si-NC with different surface 
passivations or embedded in a Si02 matrix have been car- 
ried out within Density Functional Theory (DFT) 9 . In 
the case of the embedding medium we have also anal- 
ysed the role played by disorder considering a crystalline 
phase or an amorphous glass. For the crystalline struc- 
ture we have considered a cubic cell (of side L~ 2.1 ran) 
of /3-cristobalite (BC) Si02 which is well known to give 
rise to one of the simplest Si/Si02 interface because of 
its diamond-like structure 4 . The NC/matrix crystalline 
structure has been obtained from a Si2i6 432 matrix by 
removing the 36 oxygens included in a sphere of radius 
0.6 nm placed at the center of the cubic supercell (see 
Fig.l, top panel). The result is a structure of 216 Si and 
396 O atoms with 32 Si bonded together to form a small 
crystalline skeleton with T^ local symmetry before relax- 
ation. In such core, Si atoms show a larger bond length 
3.1 A) with respect to that of the Si bulk structure (2.35 
) . No defects (dangling bonds) are present at the inter- 
face and all the O atoms at the NC surface arc single 
bonded with the Si atoms of the cluster. The dimension 
of the matrix preserves a separation of at least 1 nm be- 
tween the NCs replica, ensuring a proper shielding of the 
introduced strain . The optimized structure has been 
achieved by fully relaxing the atomic positions and the 
cell shape and size. The host matrix gets distorted after 
the relaxation (essentially due to the metastable nature 
of the BC), producing an optical behaviour similar to 
that of a completely amorphized silica, but with practi- 
cally irrelevant effects on the NC cristallinity. Therefore, 
the final system should be considered more like a crys- 
talline NC in an amorphous Si02- 

Together with the crystalline structure, the complemen- 
tary case of an amorphous silica (a-Si02) has been con- 
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FIG. 1: Stick and ball pictures 
of the final optimized structure 
of Si32 in /3-cristobalite matrix 
(top panels) or in amorphous 
glass (a-SiC>2) (bottom panels). 
Red (dark gray) spheres repre- 
sent the O atoms, cyan (gray) 
the Si of the matrix, and yel- 
low (light gray) the Si atoms 
of the NC. On the structures, 
Kohn-Sham orbitals at 10% of 
their maximum amplitude are 
reported. In white are the 
sum of all band-edges orbitals 
localized on the NC while in 
black is first the state de- 
localized on the matrix: a) 
NC: from HOMO to HOMO- 
15, matrix: HOMO-16, b) 
NC: from LUMO to LUMO+17, 
matrix: LUMO+18, c) NC: 
from HOMO to HOMO-8, ma- 
trix: HOMO-9, d) NC: from 
LUMO to LUMO+18, matrix: 
LUMO+19. 



sidered. The glass model has been generated using clas- 
sical molecular dynamics (MD) simulations of quenching 
from a melt using semi-empirical ionic potentials 8 . The 
amorphous dot structure has been obtained starting from 
a Si2ie0432 amorphous silica cell by removing the 40 oxy- 
gen atoms included in a sphere of radius 0.6 nm placed at 
the center of the cell, as shown in Fig. 1 (bottom panel) . 
Calculations of both the crystalline and amorphous struc- 
tures have been performed using the SIESTA code 10,11 
with core-corrected Troullier-Martins pseudopotentials 
and a cutoff of lbORy on density. No additional external 
pressure or stress were applied. Atomic positions and cell 
parameters have been left totally free to move. We note 
that the crystalline and the amorphous clusters present 
a nearest-neighbor distance of about 2.46 A, strained re- 
spect to the typical Si-bulk value (2.35 A), in good agree- 
ment with the outcomes of Yilmaz et al. 12 . 
In parallel to the crystalline and amorphous Si32/SiC>2 
systems, three other structures have been studied: (i) the 
pure matrix (in the same phases as used for the embedded 
NC calculations) , (ii) the isolated NCs as extracted from 
the relaxed NC-silica complexes and capped by hydro- 
gen atoms (Si32-H), and (iii) the NCs together with the 
first interface oxygens extracted as in point (ii) and then 
passivated by hydrogen atoms (Si32-OH). In the last two 
cases only the hydrogen atoms have been relaxed. The 
goal is to distinguish between the properties that depend 
only on the embedded NC from those that are instead in- 
fluenced by the presence of the matrix. The comparison 
of the results relative to different passivation regimes (H 
or OH groups) could give some insight on the role played 



by the interface region. As last point, the structures ob- 
tained as in point (ii) and (iii) have been fully relaxed 
in order to infer, through the comparison with the un- 
relaxed structures, the role played the strain induced by 
the lattice mismatch between Si and Si02- 
Electronic and optical properties of the relaxed struc- 
tures have been obtained in the framework of DFT, using 
the ESPRESSO package 13 . Calculations have been per- 
formed using norm-conserving pseudopotentials within 
the LDA approximation with a Ceperley- Alder exchange- 
correlation potential, as parametrized by Perdew-Zungcr. 
An energy cutoff of 60 Ry on the plane wave basis have 
been considered. 

Local fields (LF) contribution to the optical properties 
in the random phase approximation (RPA) has also been 
computed using the EXC code 14 . 



III. THE ELECTRONIC AND OPTICAL 
PROPERTIES 

When an interface is created between two different 
semiconductors, as in our case between Si and Si02, a 
valence and conduction band offset takes place. Depend- 
ing on the relative position of the band edges states in the 
two material one can have type I (straddling gap) or type 
II (staggered gap) heterojunction. Through the energy 
band character we have roughly evaluated the conduction 
and valence band offset for our Si32 NC embedded both 
in the crystalline and in the amorphous silica (see Fig. 2). 
Looking at energy bands localization we have found that 
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FIG. 2: Valence and conduction band offset calculated for the 
Si32 NC in BC (left panel) and in the amorphous glass (right 
panel) . 



all the band edges states from HOMO to HOMO- 15 and 
from LUMO to LUMO+17 in BC (Fig. la and lb, white 
orbitals) and from HOMO to HOMO-8 and LUMO to 
LUMO+18 in a-Si0 2 (Fig. lc and Id, white orbitals) 
are fully localized on the NC. The states HOMO-16 and 
LUMO+18 in BC (Fig. la and lb, black orbitals) and 
HOMO-9 in silica (Fig. lc, black orbitals) are the first 
states fully localized on the matrix, while LUMO+19 in 
silica (Fig. Id, black orbital) is localized on both the 
matrix and the NC. No pure states on the matrix exists 
at least for 5 eV above the gap; this a consequence of the 
amorphous nature of the matrix. Despite the uncertainty 
on the silica conduction band offset this configuration is 
typical of a type I hetero junction as reported in Fig. 2. 
The energy values reported in Fig. 2 are underestimated 
being DFT results. GW corrections are of the order of 5 
eV for silica 15 and of about 1 eV for NC 16 . 
The electronic properties just described determine the 
optical response of the system: we expect to have a 
low energy range spectrum dominated by the NC while 
the matrix control the higher energy optical transitions. 
What about the interface region? The application of the 
Bruggeman Effective Medium Approximation (BEMA) 17 
on Si-NCs is the natural choice to investigate the inter- 
play between the dielectric function of the NC, the Si02 
matrix, and the composite system. Also, the applicabil- 
ity of the BEMA constitutes a solid tool to test the level 
of interaction between the embedded NC and the em- 
bedding matrix. Initially, through the BEMA we have 
calculated the dielectric function of the composite sys- 
tem, e, by combining that of the host matrix, eh, and 
that of the NC, e c . This procedure is not well defined 
when applied to the crystalline NC because BC matrix 
loses its crystallinity after the creation of the NC, and 
therefore the BC is not representative of the final embed- 
ding medium. Therefore in all cases we should use the 
amorphous silica for a better application of the BEMA. 
The real parts of the dielectric functions are calculated 
from the corresponding DFT-RPA imaginary parts, £2, 
through the Kramers-Kronig relations. The comparison 



of the total dielectric functions obtained starting from 
the £2 of the pure (idrogenated) NC and that of the pure 
silica glass, with that of the a-Si3 2 embedded NCs, is re- 
ported in Fig. 3a. We note a clear mismatch between the 
curves, especially in the low-energy range. Actually, the 
BEMA is based on the assumptions that the embedded 
NC is spheric and not interacting with the Si02- The 
first assumption is essentially satisfied by construction, 
because we built the NC using a spherical cut-off. Then, 
the second assumption should be conclusive for the mis- 
match. 

In order to find the dielectric function of the non- 
interacting system, we adopt the alternative approach of 
the BEMA (inverse-BEMA), in which we obtain the e c of 
the NC from the difference between the DFT-calculated e 
of the composite a-Si3 2 /a-Si02 system and Eh of the pure 
silica. The result is shown in Fig. 3b. We note that the 
dielectric function of the non-interacting system matches 
that of the strained a-Si32-OH (the Si NC plus the first 
shell oxygen atoms around it as extracted by the ma- 
trix, i.e. the Si NC plus the interface region), especially 
in the 0-6 eV energy range. This strongly supports the 
idea that a true separation is established between the 
NC+interface system and the remaining silica matrix. 
This separation, due to the non-interacting behaviour of 
the parts, allows an excellent description of the absorp- 
tion at low energies. If, instead, through the BEMA, we 
compose the dielectric function of the strained a-Si32-OH 
NCs and that of the a-Si02 we obtain the result shown 
in Fig. 3c where, as expected, a good agreement with the 
DFT-calculatcd e of the embedded system is obtained. 
To underline the important role played by the strain that 
the matrix induce on the NC, in Fig. 3d, we compare the 
£2 obtained through the BEMA combining the contribu- 
tions of the relaxed a-Si32-OH NC and of the pure silica 
glass, with that of the a-Si32 embedded NC. The result 
is definitely worse with respect to Fig. 3c where strain 
was included. 

From the DFT-RPA optical spectra it is possible to intro- 
duce a new quantity, E° , which is defined as the optical 
absorption threshold when all the transitions with an in- 
tensity lower than 1% of the highest peak are neglected. 
By neglecting very low oscillator strength optical transi- 
tions we introduce a sort of "instrument resolution" that 
can connect theoretical and experimental results concern- 
ing the optical activity of a system. The obtained E° 
values for the embedded systems are reported in Table 
I. Here E° has been calculated for Si-NC containing 2, 
5, 10, and 32 atoms and compared with that of pure sil- 
ica. While the increasing of the NC size always results 
in a reduction in E°, the effect of the amorphization is 
to introduce an additional rcdshift (A) that progressively 
increase and becomes relevant for larger NCs. As a conse- 
quence the band gap tend faster to the bulk limit in the 
amorphous system with respect to the crystalline one. 
Finally, we investigated the optical properties of the Si32 
and a-Si32 embedded NCs. In accordance with Delerue 
and coworkers 18 that found, for Si-NC larger than 1.2 nm 
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FIG. 3: Comparison between 
the dielectric function of a) 
the a-Si32/a-SiC>2 and that ob- 
tained by the BEMA compos- 
ing the strained a-Si32-H NCs 
and the a-Si02, b) the strained 
a-Si32-OH and that obtained 
by the inverse-BEMA compos- 
ing the a-Si32 /a-Si02 and the a- 
Si02, c) the a-Si32/a-Si02 and 
that obtained by the BEMA 
composing the strained a-Si32- 
OH NCs and the a-Si0 2 , d) 
the a-Si32/a-Si02 and that ob- 
tained by the BEMA composing 
the relaxed a-Si32-OH NCs and 
the a-Si02- Real and imaginary 
parts are indicated by labels. 
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TABLE I: Calculated Eg) values (in eV) for different size Si nanocrystals embedded in crystalline or amorphous silica, compared 
with that of pure Si02. In the last line the difference between the crystalline and amorphous case are reported. 



(this is approximatively our cluster dimension) , a mutual 
cancellation of self-energy and Coulomb corrections, we 
neglect Many-Body effects. Local fields effects, instead, 
seem to have much more influence for Si/Si02 aggregates, 
as shown in 16 . Hence, we have investigated the optical 
properties of the S132 and a-Si32 embedded NCs including 
LF effects. The results are shown in Fig. 4. We note a 
strong similarity of the LF effects on the crystalline and 
amorphous systems with a progressive lowering of the 
spectra intensity as approaching the low energies. This 
matches well with the previous results that associates 
the low-energy absorption to the interface. Around the 
interface a strong charge discontinuity (inhomogencity) 
occurs, and this could reasonably explain the strong LF 
effects in this region. 



IV. CONCLUSIONS 

DFT calculations of the electronic and optical prop- 
erties of a crystalline and an amorphous 32-atoms Si- 
NC embedded in a Si02 matrix are described and com- 
pared with that of the same NC extracted by the matrix 
with or without the interface oxygen atoms and passi- 



vated with H. We find that the NC/matrix heterojunc- 
tion shows a type I character with the band edges states 
and consequently the low energy range optical spectrum 
is completely due to the Si NC plus the interface oxy- 
gens. Only higher in energy the response of the system 
is dominated by the matrix. This separation established 
between the NC+interface system and the silica matrix is 
also supported by the application of the Bruggeman Ef- 
fective Medium Approximation that also underlines the 
importance of the strain induced by the lattice mismatch 
between Si and Si02 when the optical properties are cal- 
culated. 

The optical absorption threshold for crystalline or amor- 
phous embedded NCs has also been evaluated as a func- 
tion of NC size finding that, beside the quantum confine- 
ment trend, the NC amorphization introduces an addi- 
tional redshift due to which the band gap tend faster to 
the bulk limit. 

Finally, the inclusion of local fields on the DFT-RPA op- 
tical absorption spectra of the crystalline and amorphous 
systems produces important redistributions of the tran- 
sition probabilities in both cases, while not affecting the 
absorption thresholds. This supports the idea that, for 
a proper description of the optical response of the inter- 
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FIG. 4: DFT-RPA calculated imaginary part of the dielectric 
function for crystalline (top) and amorphous (bottom) Si32 
embedded NCs, with and without LF contribution. 



face, inhomogcneity effects at the microscopic level (local 
fields effects) should be included. 
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